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The dwarf LSB galaxy population of the Virgo Cluster I. 
The faint-end slope of the Luminosity Function. 
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ABSTRACT 

The widely varying dwarf galaxy counts in different environments provide a strong 
challenge to standard hierarchical clustering models. The luminosity function is not 
universal, but seems to be strongly dependent upon environment. In this paper we de- 
scribe an automated procedure for detecting and measuring very low surface brightness 
(LSB) features in deep CCD data. We apply this procedure to large area CCD survey 
fields of the Virgo cluster. We show that there are many more faint (—10 > Mb > —14) 
low surface brightness galaxies than what would be predicted from extrapolation of 
the Virgo cluster catalogue luminosity function. Over our limited range of measure- 
ment the faint end slope of the luminosity function becomes a = —1.6. The luminosity 
function is flatter in the inner regions of the cluster than it is in the outer regions. 
Although these galaxies contribute a small fraction of the total stellar light of the 
cluster, they may contribute significantly to the mass in galaxies if they have large 
mass-to- light ratios similar to those recently measured for Local Group dwarf galaxies. 

Key words: dwarf galaxies - virgo cluster - luminosity function. 
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> 1 INTRODUCTION 
• I— I 

The Virgo cluster offers the best opportunity to study in 
$_( ' detail large numbers of galaxies (also dwarf galaxies) over a 
, small region of sky. It is the nearest (d~ 17 Mpc, Tikhonov 
et al., 2000) cluster with several hundreds of galaxies (~ 
1277 members, Binggeli et al, 1984) and the largest dom- 
inant structure of the Local Supercluster. It is an irregu- 
larly shaped cluster, with a high abundance of spiral galaxies 
among the bright cluster members and a large population 
(80% of the total known galaxy number) of dwarf galaxies 
(Binggeli et al., 1985). Its size is approximately 10 degrees 
across the sky, which corresponds to ~ 3 Mpc. The crossing 
time for this cluster is « O.IHq^ (Trentham et al., 2002) 
so that cluster galaxies have had plenty of opportunities 
to interact with each other. The huge mass of the cluster 
(M = 1.2 X IO^^Mq, Fouque, 2001) accelerates the member 
galaxies to very high peculiar velocities, so that it exhibits 
the highest blue-shift measured for any galaxy (IC3258, ap- 
proaching us at 1600 km/sec). As suggested from its irreg- 
ular structure, the cluster is made up of at least three sub- 
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clusters dominated by the bright eUipticals M87, M86, M49 
and thus it is probably not in a state of equilibrium: it's 
a very complex unrelaxed system where the central part is 
dynamically old, while the large surrounding region is not 
virialized, but still infalling (Bohringer, 1995). 

The first observations of the cluster date back to 
Mechain and Messier (late 18*/i century), who noticed a 
large concentration of nebulae in the northern wing of the 
Virgo constellation. The identification of the cluster as a 
self-gravitating system consisting of hundreds of galaxies 
followed soon after Bubble's 1923 discovery of Cepheids in 
M31. The first systematic investigation of the cluster was 
carried out by Shapley and Ames (1932). Since then the 
Virgo cluster has been of primary importance for extragalac- 
tic astronomy. It has been subject to many studies at opti- 
cal, radio, IR, and X-ray wavelengths aimed at addressing its 
galaxy population, evolution and gas content (for the most 
recent studies see Gavazzi et al, 2002 (optical & multiA); 
Van Driel et al., 2000 (radio) Tuffs et al., 2002 (IR); Shibata 
et al, 2001 (X-ray)). 

The most complete optical survey of the Virgo cluster 
was carried out by Binggeli et al. in 1984, using photographic 
plates and a visual detection method. Whilst the advantage 
of using photographic plates is the wide area of sky sur- 
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veyed 6° around the center of the cluster in this case), 
the major disadvantage is the low sensitivity: the survey is 
incomplete for objects with Mb < —14, missing the numer- 
ous galaxies that dominate the numbers in the Local Group, 
for example. Binggoli ct al. (1984) carried out an extensive 
systematic study of the cluster finally producing a catalogue 
of 2096 galaxies (the Virgo Cluster Catalogue, VCC). The 
papers produced by them contain photometry and morphol- 
ogy for the galaxies, analyze morphology and dynamics and 
study specific and general luminosity functions for the clus- 
ter (Binggeh et al.,1984, 1985, 1987; Sandage et al., 1984, 
1985, 1985a). 

Our main concern in this paper is with the numbers 
of faint Virgo dwarf galaxies that the subjective, visual de- 
tection method used by Binggeli et al. and others might 
have missed in less deep surveys. The observed gala:xy counts 
are quantified by a determination of the galaxy Luminosity 
Function, LF, (described by a Schcchtor function) and par- 
ticularly from the value of its faint end slope, a. After correc- 
tion for incompleteness, Binggeli et al. found a value for the 
faint end slope of q ~ —1.35. Following on from the Binggeli 
et al. survey, Impey et al. (1988) used a photographic ampli- 
fication technique to reach lower surface brightnesses. They 
found numerous additional LSB dwarf galEixies not included 
in the VCC and derived a value of a ~ —1.7. Using the same 
technique they also obtained a steepening in the faint end 
slope of the Fornax cluster from -1.3 to -1.55 (Bothun et al., 
1991). 

A more recent study of the Virgo cluster by Phillipps et al., 
1998 (again using photographic plates) produced a steeper 
slope again (q ~ —2.0). Phillipps et al. used a statistical 
method to obtain the luminosity function, subtracting faint 
galaxy number counts of fields outside of the cluster from 
those containing the cluster. This technique is very delicate 
because two large numbers are being subtracted away from 
each other to leave a small residual. There arc also largo vari- 
ations in the background counts (Valotto et al. 2001). These 
rather steep values for the faint end slope of the luminosity 
function of Virgo are consistent with some values obtained 
for other clusters. For example Bernstein et al. (1995) found 
a value of q ^ —1.4 for the more distant Coma cluster while 
Kambas et al. (2000) found a very steep slope of a ~ —2.0 
for the nearby Fornax cluster. 

These values differ, however, from recently determined 
values of the faint end slope of the field galaxy population 
derived from the extensive 2dF and Sloan surveys, where 
a ~ -1.2 (Cross et al., 2001, Blanton et al., 2001). The field 
galaxy luminosity function also corresponds very well with 
that obtained for the Local Group, a ~ —1.1 (Mateo, 1998). 

A major concern with the measurement of the dwarf 
galaxy content in different environments, is the wide range 
of data and detection methods used by different groups and 
how this affects the derived value of the faint end slope. 
We have recently looked at fields in different environments 
(Sabatini et al, JENAM 2002, Roberts et al., in prep.) using 
data in exactly the same conditions to that described below 
for Virgo. Trentham et al. (2002) found that the luminosity 
function is strongly dependent upon environment and, as 
a preliminary result, we find very few dwarf galaxies in the 
field and around isolated galaxies, if compared with numbers 
in Virgo. 

It is clear that an environmental dependency of the lu- 



minosity function places new and challenging constraints 
on current theories of galaxy formation. For example stan- 
dard Cold Dark Matter (CDM) theories predict steep faint 
end slopes in all environments (Bullock et al., 2000). Addi- 
tions to the theory, such as supcrnovac driven winds (Dekel 
and Silk, 1986) or the influence of a reheated intergalactic 
medium (Efstathiou, 1992) have been used to suppress dwarf 
galaxy formation globally. These arguments cannot be used 
to explain the large differences in the dwarf galaxy popu- 
lation content in some environments. Recently Tully et al. 
(2002) have suggested that dwarf galaxy numbers in differ- 
ent environments depend on whether the structure in which 
they reside is formed before or after rcionization ('squelch- 
ing'). They compared the Virgo cluster, which is dwarf rich, 
to the dwarf poor Ursa Major cluster. From their simulations 
they suggest that Virgo was assembled with its dwaxf galaxy 
population before reionization and Ursa Major after, so that 
dwarf galaxy formation was suppressed. There is though, as 
they say, only 'qualitative' agreement between the model 
and the observations. Moore et al. (1999) offer more of a 
nurture rather than a nature explanation to the problem. 
Their galaxy 'harassment' model provides a mechanism for 
the formation of dwarf galaxies in the cluster environment 
through disruption of discs due to tidal interactions with 
the cluster and individual galaxies. The nature of the dwarf 
galaxy population is then very dependent on the dynamical 
properties of the cluster (crossing time, galaxy number den- 
sity). The Virgo cluster has a very short crossing time (0.1 
H-^) compared to Ursa Major (0.5 H'^) (Tully et al. 2002). 
Thus the Virgo cluster galaxies have had more opportunities 
to become 'harassed'. In a study of the dynamical proper- 
ties of Virgo cluster galaxies Conselice et al. (2001) describe 
evidence supporting the idea that the dwarf galaxies arc dy- 
namically distinct from the older cluster galaxy population 
and have thus been derived from an infalling population as 
envisaged by the Moore et al. harassment model. 

An accurate derivation of the dwarf galaxy population 
as a function of environment is the only way to distinguish 

between these models. We need a detailed comparison of the 
luminosity function derived consistently for galaxies in dif- 
ferent environments with those predicted by CDM, squelch- 
ing and harassment. This is the first of a series of papers in 
which we will study in more detail the dwarf galaxy pop- 
ulation identified in Virgo and then we will derive the en- 
vironmentally dependent luminosity function of galaxies for 
comparison with the models. 

Dwarf galaxies arc extremely difficult to detect because 
they generally have both very low luminosity and surface 
brightness. To detect objects like this we need an optimum 
filter to enhance the signal and a set of selection crite- 
ria that preferentially selects cluster members rather than 
background galaxies. We also require an automated and re- 
peatable procedure that can be applied consistently to data 
taken from a large area survey. 

The paper is or- 

ganized as follows: in section 2 we present the Isaac Newton 
Telescope (INT) Wide Field Camera (WFC) survey data 
(http:/ /www.ast.cam.ac.uk^wfcsur/indcx.php). In section 
3 we discuss the methods used to determine cluster mem- 
bership and to assess possible background contamination. In 
section 4 we describe the detection algorithm (see also Saba- 
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Figure 1. Plot of the Virgo Cluster members as from the VCC. 
Overlayed are the areas covered by the INT WFCS. 

tini et al., 1999; Scaramella et al., in prep.). The luminosity 
function and first results are presented in section 5. 



2 DATA 

The data we use are part of the INT WFC survey, that is a 
multi-colour CCD based wide field survey covering an area 
of ~ 200 deg^. The camera is mounted at the prime focus 
of the 2.5m INT on La Palma, Canary Island, it consists 
of 4 thinned EEV 2kx4k CCDs with pixel size 0.33" and 
has a field size of 34.2'x34.2' (neglecting the 1' inter-chip 
spacing). The data for the Virgo survey were acquired dur- 
ing observing runs in Spring 1999 to 2002 and consist of two 
perpendicular strips of B and I band CCD images extending 
from the center of the cluster (identified as M87) outward 
for 7 and 5 degrees respectively (see fig 2). The total area 
covered is ^ 25 deg^ and the average sky noise corresponds 
to ~ 26 B mag/ sq arcsec. With the use of the techniques we 
describe in sec. 4, these deep data allow us to study the pop- 
ulation of dwarf galaxies down to very faint limiting central 
surface brightnesses (~ 2QB mag/sq axcsec) and absolute 
magnitudes {Mb ~ —10) (see below). The results presented 
in this paper make use of the East- West strip, that includes 
M87. 

The data were preprocessed and fully reduced using the 
Wide Field Survey pipeline. This includes de-biassing, bad 
pixel replacement, non-linearity correction, flatfielding, de- 
fringing (for I,Z bands) and gain correction. The photomet- 
ric calibration makes use of several (5-10 per night) standard 
stars and the zero-points are accurate to 1-2%. For the de- 
tails see http://ast.cam.ac.uk~wfcsur/pipeline.html. World 
Coordinate System (WCS) information is embedded in the 
reduced images FITS header and the astrometric calibration 
errors are < 1 arcsec. The median seeing was 1.9 arcsec. 



3 MEMBERSHIP DETERMINATION AND 
BACKGROUND SUBTRACTION. 

Membership determination is one of the crucial and con- 
troversial problems in the study of gala^Ky populations in 
nearby clusters (Valotto et al., 2001). The properties of the 
background galaxies need to be studied in detail in order to 



minimize contamination in the sample. With the aim of de- 
riving selection criteria which will enable us to separate clus- 
ter from background galaxies, we have carried out numerical 
simulations of the galaxy population in order to identify dis- 
tinguishing properties of member and non-member galaxies. 
The final goal is to find selection criteria that minimize the 
background contamination in the sample and at the same 
time maximize the number of Virgo cluster members se- 
lected. We have then tested the validity of the derived se- 
lection criteria in two ways (see ahead): 

(i) by considering the fall in surface density of our detec- 
tions with increasing distance from the cluster centre. 

(ii) by measuring the number of background galaxies de- 
tected, using our selection criteria, on fields outside the clus- 
ter. 

3.1 Numerical Simulations of a cone of universe 

The effects of rodshift on the properties of a galaxy are 
mainly the following: 

(i) surface brightness dimming due to distance (with a 
dependency of (1 -|- 0)* -|- K-corrections ); 

(ii) change in apparent size with distance. 

As a result, in principle, an intrinsically bright distant 
galaxy may appear eis a faint nearby low surface brightness 
one and could be included in the cluster members catalogue. 
For example a L, galaxy {Mb ~ —20, /io=21.7 Bfi) at z—0.2 
has an apparent total magnitude of 19.5 and a /io of 23 and 
could be mistaken for a cluster low luminosity LSB galaxy. 
Its scale size however would be small (of order 2 arc sec) 
while a cluster dwarf with the same apparent magnitude 
would have a surface brightness of 24.3 Bfj, (because of the 
surface brightness magnitude relation, see below) and hence 
a scale size of order 4 arc sec. Thus in principle cluster galax- 
ies can be distinguished by their larger sizes and fainter sur- 
face brightnesses at a given magnitude. 

In order to quantify the trend of the apparent parame- 
ters of a galaxy with redshift we have carried out some nu- 
merical simulations: we adapted a code (Esslinger-Morshidi 
Z., PhD Thesis, 1997) that populates a cone of universe and 
detects the galaxies after applying given selection criteria. 
The aim is to obtain an estimate of the completeness and 
contamination ^ of our sample for different selection criteria. 
We then use the most efficient selection criteria that prefer- 
entially select Virgo cluster galaxies. The code is composed 
of two parts: 

1) input galaxy creation; 

2) detection of the galaxies. 

The cone of universe is randomly uniformly ^ populated in 

t In what follows by completeness we mean the ratio of the num- 
ber of Virgo members selected according to the selection criteria 
to total Virgo members in the simulation. In the same way, the 
contamination is given by the ratio of the number of background 
galaxies selected to total number of galaxies (Virgo members + 
background) 

t Density fluctuations are not considered here, as the main goal 
of the simulation is to compare general apparent properties (mor- 
phology and photometry) of the background population with the 
Virgo Cluster galaxies. 
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space using a given cosmology, a Luminosity Function and a 
Surface Brightness Distribution for the galaxies. The output 
of this first part of the code is a catalogue of objects with 
given proper distance, apparent magnitude, central surface 
brightness and scale length (all objects arc assumed to have 
exponential profiles). In the second part of the code, the 
given selection criteria axe applied to the catalogue and the 
final output is a list of the objects that satisfy the selection 
criteria along with their photometric properties. 

In our simulations wc chose Ho — 75 km s^^ Mpc~^, 
a flat universe with = 0.3 and Qa = 0.7 (Lahav et al, 
2002), and we analyzed a cone from redshift 0.001 to 1.5 and 
size 10 X 10 sq degrees 

The luminosity function we adopted for the field is the 
one obtained by the 2dF Galaxy Redshift Survey (Madgwick 
et al., 2001) with the following parameters: 

• Mb* = -19.79 

• a = -1.19 

• 0* = 0.00157 h^ Mpc"^ 

and the range of input magnitudes for the galaxies is: -23 to 
-10. 

The surface brightness distribution was chosen to fol- 
low the relationship between Surface Brightness and Total 
Magnitude (for recent confirmation see Blanton et al., 2001). 
Here we use the one given in Driver, 1999: 



Me ~ [(0.67 ± 0.08)Ms - (33 ± 6)] 



(1) 



where /ie is the average surface brightness within the effec- 
tive radius and for an exponential profile: 



Me 



' ^J.o + 1.15 



(2) 



Using the same code with the same area (but keeping 
a fixed redshift of z=0.0040 ± 0.0003, corresponding to the 
distance of the Virgo cluster, Russell et al. 2000) we gen- 
erate an artificial cluster. The Luminosity function used to 

describe it is normalized to the one in Sandage et al. '85 in 
their survey of the Virgo cluster: 

• M, = -21.4 B mag 

• 0* = 2.36 Mpc"^ 

We leave the slope a as a free parameter, with the 
purpose of running simulations with different values (a = 
-1.0,-1.2,-1.4,-1.6,-1.8,-2.0). The surface brightness 
distribution chosen to describe the Virgo cluster is the Sur- 
face Brightness Magnitude relationship as derived by Impey 
et al. (1988) for Virgo: 



Ho = 0.67M + 32 ±2 



(3) 



The output catalogues produced by the code allow us to 
study the behaviour of apparent central surface brightness, 
scale length and total magnitude for simulations of the two 
different environments. Wc can then estimate what is the 
percentage of Virgo members and background gala^xies for 
each bin of scale length h and central surface brightness 
Mo. In fig. 2 we show how the percentage contamination by 
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Figure 2. Plot of background contamination as a function of 
110 (mag/sq arcsec) and h (arcsec). The numbers refer to the 
simulation where the LF has a faint end slope of -1.4 for the 
Virgo Cluster. 



background galaxies depends on the scale length and cen- 
tral surface brightness of the galaxy, for the simulation with 
a = —1.4 for the Virgo LF ^. Fig. 3 shows the histogram of 
scale lengths for all simulated galaxies. The dashed dotted 
line refers to cluster members and the filled line to back- 
ground galaxies. Both figures clearly show that the best 
discriminating property to distinguish cluster galaxies from 
background ones is the scale length: as a general trend the 
background galaxies appear to have mainly smaller sizes and 
this makes it possible to separate them from Virgo members 
(as discussed earlier). Selecting galaxies with scale lengths 
h > 3 arcsec ensures that we maximize the detection of 
cluster galaxies compared to background galaxies (from fig. 
2 contamination within this selection criteria is kept under 
50 %) . Being primarily interested in the low surface bright- 
ness population of the cluster we then restrict the selection 
to galaxies with mo > 23. 

Simulations with LF faint end slopes for the cluster 
steeper than that of the background give a better discriminar 
tion between cluster and background galaxies as the cluster 
slope increases. This is illustrated in tables 1,2,3. 



In summary our simulations have shown that in order to 
maximize the ratio of Virgo members to background galax- 
ies we should use a selection criteria of Mo > 23 and a scale 
length h > 3 arcsec. On a practical note we found that typ- 
ically the seeing on our frames was a poor 2 arc sec. Con- 
volving a 3 arc sec scale size galaxy with this seeing leads 
to a measured scale size closer to 4 arc sec. Thus in practice 
we selected objects with measured scale sizes greater than 
3 arc sec which, because of the template sizes used in the 
method described below, means a minimum scale size of 4 



§ The large aperture is necessary because otherwise the local uni- 
verse is underrepresented, due to the very small volume. Moreover 
this is also roughly the size of the cluster 



^ We show results form this simulation because this is the flattest 
LF faint end slope we would expect if we don't find any additional 
galaxy to the VCC ones 
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Figure 3. Histogram of Scale Lengtli for both simulations: the 
dashed dotted line refers to Virgo members and the filled line to 
background galaxies. 

Slope of Virgo LF = -1.0 
Scale Length cut-off Completeness Contamination 

2 93% 62% 

3 83% 31% 

4 72% 15% 

Table 1. Completeness and contaminations for different choices 
of the minimum scale length selected. 

Slope of Virgo LF = -1.4 
Scale Length cut-off Completeness Contamination 

2 88% 6% 

3 70% 2% 

4 55% 0.8% 

Table 2. Completeness and contaminations for different choices 
of the minimum scale length selected. 

Slope of Virgo LF = -2.0 
Scale Length cut-off Completeness Contamination 

2 82% 0.02% 

3 60% 0.00% 

4 44% 0.00% 

Table 3. Completeness and contaminations for different choices 
of the minimum scale length selected. 




Figure 4. At the top we show the output of SExtractor for an 
image with a LSB galaxy from our catalogue in the centre: as can 
be seen, the program classifies galaxy is identified as a group of 
small objects rather than as an object as a whole. The detection 
in this case is difficoult because of the poor pixel-to-pixcl signal 
to noise ratio for this object. As a comparison at the bottom we 
show the 3D output image from our algorithm for the same field. 
In this image the galaxy appears to be far above the noise level 
and thus easily detectable as a single object. 

chose a number of random fields to compare our predicted 
number of background detections with that of the model. 
The number of galaxies detected by our algorithm is 4 galax- 
ies per sq deg. This number is in agreement both with what 
our numerical simulations predicted (~ 6 gal per sq deg) 
and from what we measure for the background counts to be 
in the Virgo fields. It is much less than the actual number 
of galaxies ( ~ 20 gal/ sq deg) found in the Virgo cluster 
fields. All these points are discussed further below. 



3.2 Offset fields 

As a further check to validate our selection criteria, we have 
applied our algorithm (see sec. 4) to another set of INT WFC 
data using the same detector, exposure time and filter and 
covering a region of sky at about the same Galactic latitude. 
The data we used are part of the Millennium Galaxy Sur- 
vey that is a 36 arcmin wide strip going from 9''58™'28* to 
14''46™45", J2000 (Liske et al., 2002). From these data we 



4 THE DETECTION ALGORITHM 

Looking for objects with a surface flux close to the sky noise 
requires the use of image enhancement techniques in order 
to optimize their detection. Standard detection algorithms 
for connected pixels in fact often fail on these kind of ob- 
jects because of their poor pixel-to-pixel signal to noise ratio 
(S/N). The algorithm we have developed uses frequency do- 
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main techniques and mainly consists of convolutions of the 
image with matched filters: the advantage in this case is the 
use of the total flux of the galaxy to detect it, instead of the 
very low S/N pixels at its edge. As an example of this, in fig. 
4 we show the finding chart produced by SExtractor (Bcrtin 
& Arnouts, 1996) for one of our fields: a LSB galgixy from 
our catalogue (/xo = 25.3 mag/sq arcsec, h =9 arcsec) is at 
the centre of the image. As a comparison we also show the 
output of our detection algorithm on the same image: where 
SExtractor finds many distinct objects at the position of the 
LSB galaxies, our procedure clearly distinguishes it as one. 

The algorithm wo used is a modified version of the 
method described in Sabatini et al. (1999) H (for more de- 
tails see Scaramella et al., in preparation) and it is composed 

of the following steps: 

(i) background fluctuation flattening 

(ii) removal of standard astronomical objects (like stars, 
bright galaxies, cosmic rays etc.) 

(iii) image convolution with matched filters that are op- 
timized to enhance faint fuzzy structures 

(iv) candidate classification. 

LSB and dwarf galaxy candidate identification is per- 
formed on the final convolved image by means of select- 
ing all peaks that arc significantly (sec sec. 4.3) above the 
residual noise fluctuations. The whole procedure (that in- 
cludes masking, filtering, detection and measurement) is au- 
tomated and its efficiency has been tested using artificial 
galaxies added to real images (see sec. 4.5). 

Techniques involving convolution with matched filters 
have been used before (see for example Armandroff et al. 
1998 or Flint et al. 2001) but the innovations of our method 
are mainly: 

1) a generalization of the use of the convolution: we make 
use of a set of filters with several different sizes and we then 
combine all the convolved images in just one final signifi- 
cance image where objects of different scale lengths are em- 
phasized at the same time (see par 4.3). 

2) an optimization of the detection and photometry opera- 
tions, that are performed at the same step using the final 
significance image (see par 4.4 and 4.5 for more details). 

4.1 Preliminary processing 

A very important issue when optimizing the detection of 
LSB galaxies is ensuring that the sky is as fiat as possible 
across the image. Although the INT WFCS pipeline pro- 
vides fiat fielded images, wc use an automatic tool from 
the package Sextractor in order to remove possible residual 
background fluctuations. The technique consists of creating 
a map of the background sky by interpolation of the mean 
values of pixels in a grid of the original image. The grid 
size is chosen such as to preserve the biggest scale length we 
want to detect with our convolution technique. Noise reduc- 
tion gained in this way is about 6%, but a flatter background 
allows an easier and loss contaminated application of filters 
and a subsequent improved measure of source fluxes. 

II Data management and processing is performed using the appli- 
cation package language IDL (Copyright Research System, Inc.) 
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Figure 5. Sextractor detections for one of our typical images 

shown in the plane (isophotal area, surface isophotal flux weighted 
by peak flux). It is clearly possible to discriminate the stellar 
locus, a region of saturation and a region occupied by diffuse 
objects such as galaxies 

4.2 Removal of stars and standard objects 

Before convolving the image with the set of filters, we re- 
move all the 'standard' astronomical objects, such as stars, 
very bright galaxies, satellite tracks, hot pixels etc . This 
is necessary in order to minimize the contamination of the 
sample by spurious objects, i.e. objects that could simulate 
LSB galaxies when convolved with the filters. 

The masking procedure is performed in two steps in 
which different kinds of objects are considered: a first step 
is aimed at the removal of big and bright objects, such as 
saturated stars or bright galaxies, and a second one aimed 
at small stars. Although Sextractor has the option of giving 
a cleaned output image, the removal is not always effective 
and many stellar halos remain. We have thus written our 
own code for bright or saturated objects and use the Sex- 
tractor procedure for the few small stars left over after our 
cleaning. Removal is performed by masking these objects 
with the local median sky value with added Poissonian noise. 
We use SExtractor just to detect all the standard objects in 
the image and select the ones to reject using a criterium 
based on dimension (isophotal area) and peak flux (surface 
isophotal flux weighted by peak flux), that clearly discrim- 
inates a stellar locus, a possible region of saturation and a 
region occupied by diffuse objects like galaxies (fig. 5). 

Once the different regions in the plot are recognized, 
it is possible to identify the objects wo want to mask (fig. 
5) **. The stellar locus is fitted on each image so that the 
procedure is not affected by any change in the seeing value. 
Once the stellar locus is fitted we consider a line with same 
slope, but constant coefficient lower of l.bafu and mask all 
the objects that arc above this line. The object removal is 
performed masking the region with the median sky value 
plus its Poissonian noise. The size of the mask depends on 
the geometrical parameters that Sextractor gives for each 



** The theoretical trend of the stellar locus is also easily repre- 
sentable assuming a gaussian radial profile with a width equal to 
the psf 
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Figure 6. The top figure is a real image to which we added 3 ar- 
tificial galaxies of scale length (arcsec) and central surface bright- 
ness (flux/sq arcsec) as follows; (5,liTsi.y),(4,2a-gj.j^), {3,3cr^ky)- 
The 2 top left detections are real galaxies in the image. 
Bottom figure is the final output of the algorithm it is clear that 
the S/N ratio is improved and they are easily detectable. 

object (elliptical axis and Kron radius), on peak intensity 
and on the radius at which the star flux falls below la sky- 
Finally the removal of small stars left after this proce- 
dure, is carried out using the star-subtracted image given by 
SExtractor. We obtain, in this way, a 'star-cleaned' image. 

The star-masking procedure decreases the useful area 
for galaxy detection, so we take this into account when cal- 
culating numbers per sq deg in the cluster. 

4.3 The multi-scale filter 

The flrst problem when choosing the best fllter for the de- 
tection of LSB or dwarf galaxies is the determining of the 
optimal scale size for the filter - the sizes of these galaxies 
are highly variable. The necessity of selecting many different 
scales requires the use of either a very wide band pass filter 
(with the unavoidable consequence of including many dif- 
ferent kind of objects) or the application of many different 
filters where we analyze the result from each one of them, 
requiring a huge working time. Thus we decided to build 
a procedure in which we apply a combination of filters of 
different sizes but we obtain (as described below) just one 



final significance image. This final image has the property 
of having each different size emphasized in it at the same 
time. We then use it as a map of candidate positions. 

The filters have exponential profiles and are adjusted 
so that the convolution with a constant values, such as an 
empty area of an image, gives zero as output. Each filter is 
equal to zero for r = 3a and r = 6a (which means that 
it weights as positive everything inside r = 3a (typical size 
of an exponential object of scale length 3a) and subtracts 
whatever is between r = 3a and r = 6a). In this way we 
obtain a scale selective filter: everything that is smaller or 
bigger than the filter scale size is severely dimmed. 

The cleaned image is convolved with each filter and the 
output is a series of convolved images. Each of them im- 
proves the signal to noise ratio (S/N) of objects with scale 
sizes matching the filter scale. Using these convolved images 
we build a final output image whose value in each pixel is 
equal to the maximum value assumed in the stack of con- 
volved images . The main property of this final image is 
that objects corresponding to all the different sizes of filters 
are emphasized at the same time (fig 6). 

The final image is used as a map for the positions of 
candidate galaxies; both detection and photometry are then 
performed on it. It is a significance image (as the value of 
each pixel is expressed in multiples of noise of the convolved 
image correspondent to the best matching filter) and the 
detection of the candidates is done taking all peaks above a 
3(7 threshold (see sec. 4.4). On a different array we save the 
filter scale on which each maximum for each pixel is found. 
That filter scale corresponds to the best matching scale, and 
this is our estimate of the scale length of the object. 

From a computational point of view, once the filter scale 
sizes are decided, the set of filter arrays are built just once 
and then restored in the code at the convolution step. As a 
first application, we use filters with a radial symmetry, but 
the code also allows for the possibility of using an elliptical 
symmetry, which is optimized for edge on or intrinsically 
elongated objects. 

The peak value on the output image also contains all 
the information we need for the photometry of the object, 
as shown in the following section. 

4.4 Candidate classification: parameter estimation 

Given the properties of the filters (see section 4.3) it can be 
shown that when the scale size of the filter H matches the 
scale size of the galaxy h (which is the filter scale for which 
the value of the convolution is found to be maximum), the 
value of the convolution integral is: 

I = J fod{F^e-Ti)2TTrdr ^ ^TTfoh'^ (4) 

where foe~~ is the assumed exponential profile of the galaxy 
and F,e~TT is just the positive part of the filter, as the 
object has zero fiux on the negative part. Knowing the best 
matching scale, which corresponds to our best estimation of 
the object scale length, we can estimate the central original 
flux as: 

tt The intensity on each convolved image is measured as multi- 
ples of the convolved image noise, so that the measure of intensi- 
ties on different convolved images can be compared 
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fo 



(5) 



and calculate the central surface brightness. 

In the ideal case of a poissonian distribution for pixels 
in the image and assuming pixels to be uncorrolatcd, we 
can derive a simple relationship between the noise of the 
original image (at) and the noise in each convolved image 
{(To)- Within these assumptions, for a convolution with a 
simple n X n box filter, the noise reduction is: 



(To 



(6) 



The relation still holds for the exponential filters and in 
general we can then write: 

oo = k{U)o, (7) 

where the coefficient k is different for each filter and is re- 
lated to the filter scale size H. 

Using cqs. 5 and 7 we can then show how the detection 
threshold on the final convolved image (t = Ntao) relates to 
a threshold in the minimum central surface brightness which 
is detected for each scale length: 



/<), 



min — ^^,'2 p, — 



(8) 



However the assumption that pixels arc uncorrclated is 
not completely true (because of PSF) and these relationships 
need to be calibrated for our data. In order to calibrate them 
we used artificial galaxies added to real images. The results 
of these simulations are shown in the following section. 

4.5 Artificial galaxy simulations 

In order to test the efficiency of the method, we ran simula- 
tions using artificial galaxies added to real images. The use 
of artificial galaxies allows us to test the algorithm, explor- 
ing systematically the parameter space of scale length (h) 
and central surface brightness (/Jo) and thus total magni- 
tude. The added artificial galaxies have exponential profiles 
(as observations indicate that this is the best representation 
for dwarf and LSB galaxies profiles. Davies et al., 1988), are 
convolved with a gaussian that simulates the seeing and are 
added to the real data with their poissonian noise. The use 
of these simulations allows us to determine the efficiency for 
detection and photometry of the objects. 

In fig. 7 we plot the efficiency of detection as a function 
of scale length and central surface brightness of the artificial 
galax;ies, where difi^erent symbols refer to different efficien- 
cies. The efficiency is very high (90 — 100%) over almost all 
of the simulated region. As expected the efficiency drops for 
small and faint objects because of their low S/N. We now 
have an estimate of completeness and contamination of our 
detection method and we can correct for these effects. The 
objects that we can detect are just the objects we want to 
look for in Virgo. The faintest and smallest galaxies we are 
able to detect correspond to Mb ~ —9.5 and the brightest 
and biggest to Mb ~ —14.5. 

Once detected, and knowing the best matching filter 
scale, we can determine the central surface brightness (using 
eq 5) and thus the total magnitude of the objects. In figure 
8 we show the difference between input values and recovered 
ones for the total magnitudes of the simulated gala^xies. The 
efficiency in recovering the object's scale length is strongly 
dependent on how big the gap is between the scale size of 
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Scale Length 

Figure 7. Detection efficiency plotted as a function of scale 

length h and fio- The CSB (^o) is measured as a multiple of the 
sky noise and in our images this correspond to: la ^ 26.1mag/sq 
arcsec, 2cr ^ 25.3mag/sq arcscc,3(T ~ 24.9mag/sq arcscc,4CT ~ 
24.6mag/sq arcsec, 5a ~ 24.3mag/sq arcsec,10<7 ~ 23.6mag/sq 
arcsec. As shown in the plot, different symbols refer to different 
efficiencies. 
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Figure 8. The measured magnitude minus input one is plotted 

against input total magnitude. Different symtjols refer to different 
values of the /io- The biggest errors are for very faint objects. 



different filters. In order to obtain a good sampling of the 
scale lengths that we expect for dwarfs in Virgo, we decided 
to use the following 2,3,4,5,6,7 and 9 arcsec filter scales. 
Although we used smaller filters the final minimum scale 
size for objects in our Virgo sample is 4 arc sec (see sec. 3.1 
for comments on this). 

Although we do have some scatter in the recovered scale 
length and central surface brightnesses, the two compensate 
to give an estimation of the total magnitude with mean error 
of ±0.5. 

The standard way of measuring photometric parame- 
ters is by fitting to the radial surface brightness profile. The 
detection and fitting arc two distinct operations. Using our 
method we maximize our detection efficiency by detecting 
the entire image, rather than just its poor S/N edge, and we 
also obtain the best fitting parameters at the same time. 



© 2002 RAS, MNRAS 000, 1-12 



The dwarf LSB galaxy population of the Virgo Cluster I. 



9 



5 FIRST RESULTS: THE LUMINOSITY 
FUNCTION 

5.1 The radial number density profile of the 
cluster 

Clusters are very interesting regions to study in order to 
understand the role played by the environment on galaxy 
formation and evolution. They are also important because 
it might be expected that the cluster population surface 
density decreases with radius from the cluster centre and 
eventually merges into the field. This would imply that the 
properties of galaxies in the outskirts of the cluster must be 
linked to the ones of the field population. 

The sample of galaxies that wo have detected with our 
technique extends from the centre of the cluster (identified 
as M87) outward for 7 degrees. Over the resulting area of 
14 sq degrees we have identified 105 new extended dwarf 
LSB previously uncatalogued (see VCC, Impey et al. 1988, 
Trentham et al. 2002). 

Before discussing the implications on the luminosity 
function, though, we need to demonstrate that we have a 
Virgo cluster sample rather than a sample contaminated by 
background (or foreground) objects. 

Background contamination has been one of our main con- 
cerns. In fig. 9 we have plotted the surface number density of 
our detections against cluster radius. This is the raw data: 
there are no corrections for contamination or completeness 
in this plot. As expected for a cluster member population 
the density decreases when going further from the centre 
and eventually drops to an almost constant value close to 
zero, approaching the cluster edge. If our sample were highly 
contaminated by background galaxies, we would expect an 
almost flat distribution of galaxies, not one dependent on 
distance from the centre of the cluster. 

An exponential plus a constant (background) fit to the 
distribution gives a scale length of 2.2° ± 0.2° (at the dis- 
tance of the Virgo Cluster this corresponds to 0.7 Mpc) and 
a background galaxy density of 5 ± 1 gal per sq deg. This is 
consistent within the errors with the value of 4 gal/sq deg 
found for the background in the offset fields described in 
section 3.2 and can thus be considered as the non-members 
contamination in our sample. In summary, although there 
might still be some contamination by background galaxies, 
we demonstrated it is minimal (see our numerical simula- 
tions, the offset fields number counts and fig 9) 

This radial distribution can be compared with that of 
the bright galaxies. We have defined a Dwarf-to-Giant Ratio 
(DGR) which we will use in subsequent papers for compar- 
ison with different environments. This is the ratio of dwarf 
galaxies, definided as those with — 14 < Mb < —10, to giant 
galaxies with Mb < —19. We use this quantity because in 
some environments there are too few galaxies to construct 
a luminosity function. 

In fig 10 we show this ratio as a function of distance from 
M87. Interestingly this ratio remains rather flat with a me- 
dian value of ^ 20: the dwarfs number and giants one decline 
with clustercentric distance with the same scale length, re- 
sulting in a constant DGR. 

Sabatini S., Roberts S. and Davies J. (JENAM, 2002) have 
shown evidence that the DGR is about 4 for the field popu- 
lation. This is just about the number you would obtain with 
our selection criteria if observing the Milky Way from the 



Distance from M87 (deg) 

Figure 9. In this figure we plot the number density profile of 

the galaxies detected witliin our selection criteria (h> 4, > 23) 
versus distance from the cluster centre (identified as M87). As 
expected for a cluster member population the density decrease 
when going further from the centre and eventually drops to almost 
zero approaching the cluster edge. 



distance of the Virgo cluster. Out of the group of galaxies 
that Mateo (1998) assigns to the Milky Way just the dwarfs 
Sextans Fornax and Sagittarius would meet our selection 
criteria. Again this would give a value for the DGR of 3, 
lower than in Virgo. As most dwarf galaxies in clusters are 
probably not bound to individual giant galaxies we can also 
do the same test for all the dwarfs in the Local Group, in 
order to compare it with the Virgo Cluster. Again we obtain 
a DGR of 4. 

Also, interestingly, the Virgo cluster DGR does not appear 
to smoothly blend into the field; if it did we would expect 
DGR to gradually decrease to the value in the field. Note 
also that we must be close to the 'edge' of the cluster be- 
cause the galaxy counts are about the same as those in our 
offset fields. Thus the Virgo cluster environment seems to 
be very different to that of the field even in its most outer 
regions. 



5.2 The faint end slope of the Luminosity 
Function. 

In this paper we are primarily interested in the number of 
dwarf galaxies in the cluster. When sufficient galaxies are 
available this is found by fitting to the faint end of the lu- 
minosity function. Our method enables us to detect galaxies 
with the following range of absolute magnitudes: Ms =-10 
(h™„ = 4", Atr" = 26 Bii) to Mb = -15 (p^"" = 23, 

hmax — 9). 

An important check is to compare our derived magni- 
tudes with galaxies common to previous catalogues of the 
Virgo Cluster. We have 143 galaxies in our sample that are 
listed in Trentham et al. 2002 and we show a plot of our 
measured apparent magnitudes against the Trentham et al. 
measurements in fig 11: a linear fit to this plot results in a 
slope of 1.02 and a constant value of -0.22. Our magnitudes 
tend to bo slightly brighter then Trentham's, but the differ- 
ence lies within the errors we expect in our measurements 
(see fig 8). Galaxies listed in Trentham et al. 2002 that are 
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Figure 10. Dwarf to giant ratio as a function of distance from 

M87. Here wc define dwarf tlie galaxies witti — 14 < Mb < — 10 
and giant tiie ones with Mg < —19. Tiie pealc at ~ 5° is mainly 
due to a dip in the distribution of giant and not to an excess of 
dwarfs. 



Figure 12. In this figure we plot the faint end of the luminosity 

function. Open circles arc counts for the raw data, while filled 
points arc corrected for detection efficiency and background con- 
tamination. Squares are points from the VCC luminosity function 
with slope a = —1.35, renormalized for our survey area. 




Figure 11. Comparison between apparent magnitudes as mea- 
sured with our method (x axis) and those in Trentham et al. 
2002. A linear fit to the data gives a slope of 1.02 and a constant 
of -0.22. 



not in our sample have been checked and they are either 
too bright to fall within our magnitude range or they lie in 
masked regions of the images (in our calculations we take 
into account the area lost due to the removal of stars). 

Our derived luminosity function is shown in fig. 12. We 
show the raw data and the corrected data. The number of 
objects detected in each bin of h and fj,o has been corrected 
for: 

(i) background contamination (as obtained from our nu- 
merical simulations of a cone of the universe). 

(ii) incompleteness (as estimated by the detection effi- 
ciency of the algorithm with artificial galaxies). 

The corrections make little difference to the numbers de- 
tected in each bin. Assuming that the drop-off in numbers 
beyond Mb = —10.5 is due to incompleteness we have fit- 
ted the luminosity function in the magnitude range -14.5 to 
-10.5. This gives a value for the faint end slope of —1.7 ±0.2 



Central Surface Brightnt 



Figure 13. Histogram of Central Surface Brightness distribution 
for our detections. 



for the raw counts and —1.6 ± 0.1 for the counts corrected 

for incompleteness. 

It is difficult to combine galaxy counts from samples 
selected in different ways. For example in Kambas et al. 
(2000) we show how the selection of three different samples 
of galaxies leads to a disjoint surface brightness distribution 
because in each case the selection criteria preferentially picks 
galaxies of a given surface brightness. The surface brightness 
distribution of this sample (fig. 13) is also peaked. 

Even so in fig. 12 wc have also shown, for compari- 
son, the number counts for the VCC tt. Including the VCC 
data and fitting the faint end slope between -16 and -10.5 
we obtain a value of -1.8 for a The values that we ob- 



M These can also be used to compare our data with the number 
counts obtained in Trentham et al 2002, as his findings fit very 
well with the Luminosity Function of the VCC (see their fig 5). 

We can't write an error on this value because we don't have a 
completeness function for the VCC catalogue. 
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Absolute B magnitude 

Figure 14. In this figure we plot the faint end part of the lu- 
minosity function separating the cluster into an inner (distance 
from M87 < 0.8°) and outer (distance from M87 > 1.6°) region. 
Dashed-dotted line refers to the former, while filled line to the 
latter. Although the numbers give a poor statistics, a K-S test 
shows that the two distributions are different with a 90% proba- 
bility. A fit to the same range of magnitudes as for the total LF, a 
steeper value for the faint end slope of the LF in the outer region 
compared with the inner one (—1.8 ± 0.2 and —1.4 ± 0.2 ) 



tain are in agreement with a steepening of the LF when in- 
cluding in the sample the contribution from faint and LSB 
galaxies. The original VCC value of -1.35 (Sandage et al., 
1984) had already been brought into question by Binggeli et 
al. (1988) who suggested a steepening to -1.7 when includ- 
ing a possible dwarf population that had escaped detection. 
The final slope we obtain is not as steep as the one ob- 
tained by Phillipps ot al. (1998) using the background field 
subtraction method. However, in order to compare our re- 
sults with Phillipps et al., we should consider our raw data 
counts (as they didn't make any correction for incomplete- 
ness) that result in a slope of — 1.7±0.2. Neglecting the last 
point in their LF, that might bo highly background contam- 
inated (S. Phillipps, private communication), their slope is 
~ —1.9 ± 0.15. This then is consistent with our result. 

We can also calculate a separate luminosity function 
for the inner and outer region of the cluster (see fig. 14). 
Dashed-dotted line refers to the former, while the filled line 
to the latter. We have chosen these regions because they 
correspond roughly to that part of the cluster that is within 
the virial radius and that which is outside it. Galaxies within 
the virial radius should have been much more affected by 
interactions with other cluster galaxies than those in the 
dynamically unrelated outer regions (Bohringer, 1995). A fit 
to the same range of magnitudes as for the total luminosity 
function results in a steeper value for the faint end slope in 
the outer region compared with the inner ( — 1.8 ± 0.2 and 
— 1.4 ± 0.2 ). A K-S test shows that the two distributions 
are different with a 90% probability. If confirmed by a larger 
statistics, this result is consistent with the idea that the 
faintest galaxies arc more abundant in the outer regions of 
clusters, while in the denser inner regions they have partly 
been accreted by larger galaxies or have dimmed or even 
been disrupted by tidal interactions. 



5.3 Total Luminosity and Mass 

The total light (corrected for the area of the cluster sampled 
compared to the VCC) due to the population of dwarf LSB 
galaxies of our sample is 

Lowarf = 4 X lO^L© (9) 

which is just 1/50 of that due to VCC galaxies. So this 
population of galaxies contributes only a small fraction of 
the light contributed by the bright galaxies. 

A surface brightness level of about 28 B mag/sq arcsec 
has recently been predicted for intra-cluster light from stars 
associated with intra-cluster planetary nebulae (Arnaboldi 
et al. 2002). The average integrated surface brightness from 
the galaxies we detected is fitot ~ 32 mag/sq arcsec in the 
inner region of the cluster (within 0.8° from M87). This is 
far too faint to have been previously detectable as a surface 
brightness enhancement. This, combined with the planetary 
nebulae data, may indicate that there are more even lower 
surface brightness structures to be discovered between the 
galaxies. 

Despite the small contribution to the light, dwarf galax- 
ies may contribute a larger fraction to the mass. Recently 
very large mass-to-light ratios have been found for Local 
Group galaxies (see for example Mateo 1998 and particu- 
larly Kleyna et al. 2002). Assuming a mass-to- light relation 
as in Davies et al. (2002): 

M ^ ^q3.5^-0.2B (^q) 
J-i 

the total mass due to the dwarf LSB population over an area 
of ~ 14 sq deg is 

Mowarf = 1.9 X 10" Mq (11) 

Rcscaled to the total area of the cluster, this is 1/10 
of the total mass from the VCC galaxies. Given that other 
low surface brightness material associated with intra-cluster 
stars exists, it is possible that there is as much mass in the 
low surface brightness component of the Virgo cluster as 
there is in the brighter galaxies. 



6 CONCLUSIONS 

Observations of the relative numbers of dwarf galaxies in dif- 
ferent environments present a strong challenge to galaxy for- 
mation models that predict large numbers of dwarf galaxies 
in all environments. They also present a challenge to those 
models that predict global suppression of dwarf galaxy for- 
mation. The Virgo cluster is a very different environment 
from that of the Local Group and from that of the general 
field. It has a very large number of dwarf galaxies compared 
to the giant galaxy population. 

In this paper we have described a new automated tech- 
nique for finding low surface brightness galaxies on wide 
field CCD data. We have carried out simulations to ensure 
that our detection method and selection criteria enable us 
to preferentially select cluster galaxies, rather than those 
in the background. From the decrease in surface number 
density with clustercentric distance we believe that we have 
achieved this goal. 
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Over the resulting area of ~ 14 sq degrees that we anal- 
ysed, we have identified 105 new extended dwarf LSB previ- 
ously uncatalogued (see VCC, Impey et al. 1988, Trentham 
et al. 2002). The resulting luminosity function is consider- 
ably steeper than that inferred from an extrapolation of the 
data in the VCC. The cluster luminosity function appears 
to be steeper in the outer parts of the cluster than in the 
inner part, though the dwarf- to-giant ratio remains almost 
constant. Although these galaxies contribute only a small 
fraction of the luminosity of the cluster they may contribute 
significantly to the galactic mass of the cluster, given re- 
cently measured large mass-to-light for Local Group dSph 
galaxies. 

This paper has been produced using the Royal Astronomical 
Society /Blackwell Science M^^X style file. 
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